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Energy Metric for Detection 


Inputs: Luminance image = L(x,y) 

Pixel area = dx dy in deg 2 , 

Duration = dt in sec 

Compute visible contrast image = Cv(x,y) 
Visible Energy Metric : 

Ev = dx dy dt £ x y Cv(x,y) 2 deg 2 sec 

dBV = 10 log 10 ( Ev / 1 0 _ 6 ) 

Modelfest average threshold = 7 ±2 dBV 



Visible Contrast Image 

Optic Blur : Lo(x,y) = L(x,y) * 0(x,y) 

Background Luminance : 

Lb(x,y) = a(dt) (Lo(x,y) * B(x,y)) + (1 - a(dt)) B 


Contrast : C(x,y) = Lo(x.y) - Lb(x.v) 

Lb(x,y) 


Eccentricity Sensitivity : 

Cv(x,y) = C(x,y) S(x,y) 



Parameters 

Optic Blur : F(0(x,y)) = exp(- f / f 0 ), 

f = V(f x 2 + f y 2), f 0 = 12 cpd 

Background Luminance : 

F(B(x,y)) = exp(- (f / 2 ), f 1 = 2 cpd 

a(dt) = exp(- dt / 1 0 ), t 0 = 0.4 sec 

Eccentricity Sensitivity : 

S(x,y) = 1/(1 + g (1 - exp(- r / r 0 ))) , 

r = V(x 2 + y 2 ), r 0 = 5.7 deg, 

g = 4.1, 1 / (I + g) = 0.2 



Metric-Validating Model 

Visibility Image: Cv(x,y) 

Additive White Noise with 2-sided power spectral 
density 

N = a 2 dx dy dt , 

Each pixel is independently distributed as 
Normal with mean 0 and standard deviation o 

Ideal Observer detects presence or absence of 
signal in a two interval forced choice (21 FC) 
experiment. 



21 FC Model Performance 


Visibility Image: Cv(x,y) with visible contrast 
energy Ev and noise spectral density N 

Distance between observer output distributions 
divided by their common standard deviation is 

d’ = V(2 Ev / N) 

Prob(Correct) = Pc = Fz(d') - 0.5 
Estimated N = 2 Ev / d' 2 
If Pc = 0.84, d' = 1, N = 2 Ev 
Modelfest : 10 log 10 (N) + 60 = 10 ± 2 dB 



Discrimination Model 


Visibility Images: Cv(x,y,j), j = 1,M 

Additive White Noise with power spectral density 

N = a 2 dx dy dt 

Ideal Observer responds k if image j is presented 
and image k has the smallest squared distance 
d(k) to the noisy image 

d(k) = || Cv0 + N-Cv(k) || 2 

d(k) = || Cv0 + N|| 2 + || Cv(k) || 2 

- 2 (Cv(j) • Cv(k) + N • Cv(k)) 



Discrimination Model Metric 

Orthogonal Images: Cv(j) • Cv(k) = 0, j ^ k 

Same energy: Ev(j) = Ev 

Let d’ = V( Ev / N) 

Pc = J F(x)k-i f(x-d') dx , 

where F and f are the cumulative and density 
distribution functions of the standard normal. 

Also Ev = dx dy dt ^ ||Cv(j) - C|| 2 /(M-1 ) 
where C = £j Cv(j) / M 



Discrimination Model Performance 



Example: Landolt C 
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Pedestal invariance of ideal observer allows the 
orthogonal stimulus model. 



Example: Tumbling E's 



Model simulation for n = 10000 trials, d' = 1. 
95% confidence interval for Pc = 0.538 ± 0.010 
Metric prediction Pc = 0.552 




Method Considerations 


When pattern energies are similar, varying the 
contrast adds little or no uncertainty; varying 
size or blur contributes significant uncertainty. 

Practice of computing thresholds by averaging 
reversal endpoints has problems 

1) Pc at threshold is not actually known 

2) No estimate of the slope at threshold is 
provided 

3) Valuable data is effectively discarded 



Summary 


Detection metric: 

Visible contrast energy 
Approximate Discrimination metric: 
Average (M-1) squared distance from 
each visible contrast pattern to the 
mean visible contrast pattern 
Model simulation is fast 



Tumbling E Model Matlab Code 


c = s'*s ; % 4x25 times 25x4 
[u f x, v] = svd(c) ; 

f = u*(x. A 0.5) ; 

sn = 

ones (n,4)*c(l,l:4) +randn ( n , 4 ) *f 
Pc=mean ( 

sn ( 1 : n , 1 ) >max (sn(l:n,2:4) ' ) ' 

) ; 



Watson &Ahumada (2005) 

Metric Elements 



Figure 4. Elements of the component model. 


